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Abstract 



We modify the mesonic wavefunction by using a short distance scale ro in analogy with 
Hydrogen atom and estimate the values of masses and decay constants of the open flavour charm 
mesons D, D s and B c within the framework of a QCD Potential Model. We also calculate 
leptonic decay widths of these mesons to study branching ratios and life time. The results are 
in good agreement with experimental and other theoretical values. 
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1 Introduction 

Recently, we have reported a regularisation procedure to avoid the singularity by introducing a flavour 
^ \ dependent short distance scale at the origin to study the oscillation frequency of B and B mesons 
^ ■ in a QCD Potential Modelfl]. The purpose of this letter is to use the Potential Model to calculate 
the masses and decay constants of open flavour charmed mesons D, D s and B c and then to find the 
decay width and branching ratio of the same. 

If the CKM element is well known from other measurements, then f p can be well measured. If, on 
the other hand, the CKM element is less well or poorly measured, having theoretical input on fp 
can allow a determination of the CKM element. These decay constants can be accessed both exper- 
imentally and through lattice Quantum Chromodynamics (1QCD) simulations. While for f n , Jk , 
fn , experimental measurements agree well with lattice QCD calculations, a discrepancy is seen for 
the value of fn s '■ The 2008 PDG average for fo s is 273 ± 10Me^[2], about 3cr larger than the most 
precise N f = 2 + 1 1QCD result from the HPQCD/UKQCD collaboration [3], 241 ± 3 MeV. On the 
other hand, experiments and 1QCD calculations agree very well with each other on the value of fp, , 
f D (expt) = 205.8 ±8.9 MeV and f D {lQCD) = 207 ±4 MeV. The discrepancy concerning f Ds is quite 
puzzling because whatever systematic errors have affected the 1QCD calculation of fp, , they should 
also be expected for the calculation of fp s [4]. 

However, the discrepancy is reduced to 2.4a with the new (updated) data from CLEO [5,6] and Babar 
[7], together with the Belle measurement [8] and the latest PDG average is fp> s = 257.5 ± 6.1 MeV 
[9]. Lately the HPQCD collaboration has also updated its study of the D s decay constant [10]. By 



including additional results at smaller lattice spacing along with improved determinations of the lat- 
tice spacing and improved tuning of the charm and strange quark masses, a new value for the D s 
decay constant has been reported as fn a = 248.0 ±2.5 MeV. A better understanding about the decay 
properties of D and D s in this regard is quite important. 

The properties of the B c meson are of special interest [11] , since it is the only heavy meson consisting 
of two heavy quarks with different flavours. A pecularity of B c decays is that both the quarks may 
involve in its weak decays. From experimental point of view, study of weak decays of B c meson is quite 
important for the determination of CKM elements. More detailed information about its decay prop- 
erties are expected in near future at LHC and other experiments. In this letter the decay constants 
of D, D s and B c are computed and are used to study the decay width and decay time of the same. 
Rest of the paper is organised as follows: In section 2, we discuss the formalism, while in section 3 
we summarise the calculations and results. Section 4 contains summary and conclusion. 



2 Formalism 

2.1 Wave function in the model 

The wave function in a specific potential Model with its relativistic corection procedure is [TJ IT2"] : 

Ad+conf (r ) = -p^e ( C - 



where 



7r<2n 



GO, 
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N = j (2) 

'2 2 T (3 - 2e) C' 2 - lfiba 3 T (5 - 2e) C> + ^ 2 b 2 a 6 T (7 - 2e) 



C = 1 + cA^Tial (3) 
rrii + rrij 



a = ( -fia s ) (5) 



and the QCD potential is taken as 



c =l_yi-^-a.J (6) 
r' = r + r (7) 
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V (r) = -— a s + br + c 
3r 



Table 1: Values of cut-off parameter tq for different mesons 



mesons 


reduced mass 


values of Tq in Gev 1 


D(cu/cd) 


0.276f 


0.0073 


D s (cs) 


0.3682 


0.0055 


B c (bc) 


1.1838 


3.872 x 10" 1U 



2.2 The mesonic wave function at the origin 

The term ( I ^ a ) in eq.l is used to incorporate relativistic effect in a free Dirac way [12,16]. The cut 
off parameter is used to regularise the wavefunction at the origin such that the wavefunction remain 
approximately constant from ro down to r = 0. In analogy to QED, we use the cut off parameter to 
regularise the wavefunction as[l] 

r ~a e" 7 . (9) 

Unlike QED, it is flavour dependent depending on the flavours of the quark masses. 

The input parameters in the numerical calculations are used as m u = 0.336GeV, m s = 0A83GeV , 
m c = 1.550GW and m b = 4.950GW, b = 0.183GeV 2 ,c = -0.5GeV, a s = 0.39 and a s = 0.22 at 
charm and bottom mass scale. The value of the undetermined factor Ao, which appears in the series 
solution[12] is so chosen that the consistancy of the value of cAq = lGeV 3 ^ 2 is sustained with our pre- 
vius works[l,13]. With these values, we compute the short distance scale r for open flavour charmed 
mesons D, D s and B c and is shown in table 1. 

It shows as the masses of the mesons increase, r is decreased. For heavy-heavy mesons like B c it 
becomes as tiny as 10 _10 GeV A_1 which presumably infers that the relativistic effect is too small for 
this heavy meson. 



3 Calculation and results 

3.1 masses and decay constants of open flavoured charm mesons 

The decay constant with relativistic correction can be expressed through the meson wave function 
$p(p) in the momentum space [TTJEHjas 



12 f d 3 p ( E q + m q \ 1 ' 2 ( E q + m q \ 1/2 ( ^ , x p 2 
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with \ P = — 1 . In the nonrelativistic limit << 1.0, this expression reduces to the well known 
relation between f P and the ground state wave function at the origin if) (0) the Van-Royen-Weisskopf 
formula |19|. 



Table 2: The masses of heavy-light mesons in GeV. 



mesons 


masses 


experimental values 


D(cu/cd) 


1.860 


1.869 ±0.00160 


D(cs) 


1.959 


1.968 ±0.0033[9J 


B c (bc) 


6.507 


6.277 ±0.006[9J 



Table 3: The decay constants of heavy-light mesons with and without QCD correction in GeV 



meson 


our work 


other work 






fp fpc 






D(cu/cd) 


0.240 


0.225 


0.205 ± 0.085 ± 0.025 [21] [Exp] 
0.206 [22] 
0.234 [23] 
0.208 [24] 
0.201 [25] 
0.235 [26] 




D(cs) 


0.291 


0.266 


0.254 ±0.059 [21] [Exp] 
0.245 [22] 
0.268 [23] 
0.256 [24] 
0.249 [25] 
0.266 [26] 




B c (bc) 


0.435 


0.413 


0.433 [27] [Theory] 
0.470[28] 





f p = decay constants without QCD correction 
f pc = decay constant with QCD correction 



with QCD correction factor the decay constant can be written as[20] 



f P = \ ^ (0) \ 2 C\ whereC 2 = 1 - a s /n[2 - mq "^ fo^] (12) 

V M P m Q + m Q m Q 

where M p is the pseudoscalar meson mass in the ground state and can be obtained as [28] 

M p = ttlq + m-Q + (H) (13) 

where 

(H) = (£) + (V (r)) (14) 
In table. 2 and 3 we record the prediction of the model for masses and decay constants. 



Table 4: The leptonic branching ratio of D anc 


. D s mesons. 


mesons 


BRr x KT 3 


BR^ x 10~ 4 


BR e x 10~ 8 


D(cu/cd) 
Expt.[9] 
B. Patel etal.,[29] 


0.78(0.68) 
< 1.2 
0.9 


4.4(3.89) 
3.82 
6.6 


1.05(0.92) 
1.5 




BRr x 10~ 2 


BR^ x 10~ 3 


BR e x 10~ 7 


D(cs) 
Expt.[9] 
B.Patel etal.,[29] 


6.4(5.3) 
5.6 ±0.4 
8.4 


7.0(5.9) 
5.8 ±0.4 
7.7 


1.67(1.40) 
1.8 



*Values within the bracket represent the branching ratio for f p with QCD correction. 

3.2 Leptonic decay width and Branching ratio of D and D s mesons 

Charged mesons formed from a quark and anti-quark can decay to a charged lepton pair when these 
objects annihilate via a virtual W ± boson. Quark-antiquark annihilations via a virtual W^ + (W^~) 
to the l + h>(l~i> final states occur for the tt ± , K ± , DfandB ± mesons. There are several reasons for 
studying the purely leptonic decays of charged mesons. Such processes are rare but they have clear 
experimental signatures due to the presence of a highly energetic lepton in the final state. The 
theoretical predictions are very clean due to the absence of hadrons in the final state. The total 
leptonic decay width of D, D s mesons are given by 

F (flf -> |+„) = (l - ^) 2 M D „ q = d,s (15) 

These transitions are helicity suppressed; i.e., the amplitude is proportional to mi, the mass of the 
lepton 1, in complete analogy to 7r — >■ 

The leptonic widths of the charged D and D s mesons are obtained using equation. 12, employing the 
predicted values of the pseudoscalar decay constants f D and f s along with the masses of the M D and 
M Ds from our work. The leptonic widths for seperate lepton channel by the choice of m; =T)At e are 
computed. Branching ratio of D q mesons are calculated by using the relation 

BR = Txr (16) 

The life time of these mesons are To = 1.04ps and td s = 0.5ps are taken from the world average value 
reported by Particle data group (PDG-2010)[9]. The present results as tabulated in table. 4 are in 
accordance with the available experimental values. 

3.3 Weak decay of meson 

Adopting the spectator model for the charm beauty mesons system [30] , the total decay width of B£ 
meson can be approximated as the sum of the widths of 6-quark decay keeping c-quark as spectator, 
the c-quark decay with b as spectator, and the annihilation channel B+ — > l + U[(cs, us), I = e, fi,r with 
no interference assumed between them. Acordingly, the total width is written as [30] 



r (b c ->■ x) = r {b ->■ x) + r ( c ->■ x) + r (anni) 



(17) 



Table 5: Decay width(in 10 4 eV) and life time of B c meson 



meson 


T (Anni) 


r (B c X) 


r (inps) 


B c b 


1.17(1.06) 


19.17(19.06) 


0.344(0.346) 


[9] 






0.453 ± 0.041 


[30] 


1.40 


14.00 


0.47 


[31] 


0.67 


8.8 


0.75 



Neglecting the quark binding effects, the b and c inclusive widths in the spectator approximation 
are[30] 

rM= «2M (18) 

r(c^) = 5G ^y (19) 

Here we have used the model quark masses and the CKM matrix elements |V^ S | = 0.957, \V c b\ = 0.039 
from the particle data group. 

Employing the computed mass and pseudoscaler decay constant from the present study, the width of 
the annihilation channel is computed using the expression given by [28,30], 



G 2 F \V b c \ 2 P B M B ( ^ 

8tt q I Mi 1 ^ q 



T(Anm)= ^— o ^ ^ ( 1 - ^ ) C„ (20) 



Where C q = 1 for ther^ T channel and C q = 3|V^ S | 2 for cs channel, and m q is the mass of the heaviest 
fermions.The computed results of the annihilation decay rate and total decay rate are tabulated in 
table. 5. Our prediction for life time with these results lie well within the experimental value and is 
placed in table. 5. 



4 Conclusion 

In the present work, we have computed the masses of heavy-light mesons in a specific potential model 
and use a short distance scale to regularise the wave function near the origin to compute the decay 
constants. The short distance scale as well as the decay parameters are very sensitive to the strong 
coupling constant a s . It is to be noted that strong coupling constant as is used in ref.[l] cannot be 
used for D, D s and B c mesons since the results oversoots the experimental and theoretical values and 
hence the model is excluded to treat the D and B mesons together. However we use the same a s as 
is used in our previous work[13,14,15] and with a different renormalisation mass scale one can regain 
these strong coupling constants within the formula of ref. [1]. 

The analysis with QCD correction is found to be closer to experiments and other theoretical re- 
sults. In our calculation we have found ^-=1.18 (with QCD correction) and ^-=1.21 (without QCD 
correction), which is in accordance with the lattest QCD Sum rule result 1.193 ± 0.025 ± 0.007 [22] . 
PDG average ^ = 1.25 ± 0.06±[9], as well as with the recent lattice results f D Jf D = 1.164 ± 0.011 
[3,10] and f D Jf D = 1.20 ± 0.02 [32]. 

Present study on the leptonic decay branching ratio of D and D s mesons with QCD correction for 



r and /i leptonic channels presented in Table. 4 are as per the available experimental limits. Large 
experimental uncertainly in the electron channel make it difficult for any reasonable conclusion. The 
computed result within the framework of the QCD potential model for annihilation decay width as 
well as the lifetime of B c meson is also found to be well with the available datas as presnted in the 
Table. 5. Probaly, in future high luminosity better statistics and high confidence level data sets will 
be able to provide more light on the spectroscopy and decay properties of these open charm mesons. 
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5 Appendix 

Equations of calculation for kinetic energy, coulomb and confining terms of the potential to obtain 
masses of the pseudoscalar mesons. 



(H) = <^> + (V(r)> 



(21) 



'2(1' 2/ (2g) 1 -' ( ' (2g) 3 - 2 ' ( ' Pj) 2 ^ ^ ' 

_4, gftlC ,( l+ l) r (4 _ 2e) _ 2, C ,(^ +t ,(2- t )(3- ( )) r (3 _ 2£) + 



(29)— (2 9 ) 

2^,(^ 1 ,y r )(3-e)) r (5 _ 2£) + 4^ (4 _ 2£) 

-|^r ,6 - 2e) + ^filr (5 - 2e) - *££r ,7 - 2e)) 
(2c/) (2#) (20) 

(22) 

—I ' ' ••'-(! 

(23) 



(— ^> = -a c 2vr^^-^r (2 - 2e) - c^^-^r (6 - 2c) + c^Tr/^^-^r ( I 2r) 



(6r> = ^vr^^i^r (4 - 2c) + W^^g^r (8 - 2c) - Mt^^-^ (6 - 2c) 

(24) 

wherehi = hk\ 

fj,ba 
h = 



k = 



N' 



a 



(™o) 5 

d = C% 
1 

9 = — 

a 

4a s /3 



